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ABSTRACT 

A GILE is a small gamma-ray astronomy satellite mission of the Italian Space Agency dedicated to 
high-energy astrophysics launched in 2007 April. Its ^1 absolute time tagging capability coupled 
with a good sensitivity in the 30 MeV-30 GeV range, with simultaneous X-ray monitoring in the 
18-60 keV band, makes it perfectly suited for the study of gamma-ray pulsars following up on the 
CGi?0/EGRET heritage. In this paper we present the first AGILE timing results on the known 
gamma-ray pulsars Vela, Crab, Geminga and B 1706-44. The data were collected from 2007 July 
to 2008 April, exploiting the mission Science Verification Phase, the Instrument Timing Calibration 
and the early Observing Pointing Program. Thanks to its large field of view, AGILE collected a 
large number of gamma-ray photons from these pulsars (^10,000 pulsed counts for Vela) in only 
few months of observations. The coupling of AGILE timing capabilities, simultaneous radio/X-ray 
monitoring and new tools aimed at precise photon phasing, exploiting also timing noise correction, 
unveiled new interesting features at sub-millisecond level in the pulsars' high-energy light-curves. 
Subject headings: stars: neutron - pulsars: general - pulsars: individual (Vela, Crab, Geminga, 
PSR B 1706-44) - gamma rays: observations 
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1. INTRODUCTION 

Among the ~1,800 known rotation-powered pul- 
sars, observed mainly in the radio band, seven ob- 
jects have been identified as gamma-ray emitters, 
namely Vela (B 0833-45), Crab (B 0531-^21), Geminga 
(J 0633-^1746), B 1706-44, B 1509-58, B 1055-52 and 
B195H-32 (Thompson 2004). In addition, B 1046-58 
(Kaspi et al. 2000), B 0656-^14 (Ramanamurthy et al. 
1996) and J 0218+4232 (Kuiper et al. 2000, 2002) were 
reported with lower confidence (probability of periodic 
signal occurring by chance in gamma-rays of ^10"**). In 
spite of the paucity of pulsar identifications, gamma-ray 
observations are a valuable tool for studying particle ac- 
celeration sites and emission mechanisms in the magne- 
tospheres of spin-powered pulsars. 

So far, spin-powered pulsars were the only class of 
Galactic sources firmly identified by CGi? 0/EGRET and 
presumably some of the unidentified gamma-ray sources 
will turn out to be associated to young and energetic 
radio pulsars discovered in recent radio surveys (Manch- 
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ester et al. 2001; Kramer et al. 2003). In fact, several 
unidentified gamma-ray sources have characteristics sim- 
ilar to those of the known gamma-ray pulsars (hard spec- 
trum with high-energy cutoff, no variability, possible X- 
ray counterparts with thermal/non-thermal component, 
no prominent optical counterpart), but they lack a radio 
counterpart as well as a supernova remnant and/or pul- 
sar wind nebula association. Radio quiet, Geminga-like 
objects have been invoked by several authors (Romani & 
Yadigaroglu 1995; Yadigaroglu & Romani 1995; Harding 
et al. 2007) but without evidence of pulsation in gamma- 
rays, no identification has been confirmed. 

Apart from the Crab and B 1509-58, whose luminosi- 
ties peak in the 100 keV and about 30 MeV range re- 
spectively, the energy flux of the remaining gamma-ray 
pulsars is dominated by the emission above 10 MeV with 
a spectral break in the GeV range. 

AGILE (Astro-rivelatore Gamma ad Immagini LEg- 
gero) is a small scientific mission of the Italian Space 
Agency dedicated to high-energy astrophysics (Tavani 
et al. 2006, Tavani et al. 2008, in preparation) launched 
on 2007 April 23. Its sensitivity in the 30 MeV 30 GeV 
range, with simultaneous X-ray imaging in the 18-60 keV 
band, makes it perfectly suited for the study of gamma- 
ray pulsars. Despite its small dimensions and weight 
(~100 kg), the new silicon detector technology employed 
for the AGILE instruments yields overall performances 
as good as, or better than, that of previous bigger instru- 
ments. High-energy photons are converted into e+/e^ 
pairs in the Gamma-Ray Imaging Detector (GRID) a 
Silicon- Tungsten tracker (Prest et al. 2003; Barbiellini 
et al. 2001), allowing for an efficient photon collection, 
with an effective area of ~500 cm^, and for an accurate 
arrival direction reconstruction (~0.5° at 1 GeV) over a 
very large field of view, covering about 1 /5 of the sky in 
a single pointing. The Cesium-Iodide mini- calorimeter 
(Labanti et al. 2006) is used in conjimction with the 
tracker for photon energy reconstruction while support- 
ing the anti-coincidence shield in the particle background 
rejection task (Perotti et al. 2006). The AGILE /GRID is 
characterized by the smallest dead-time ever obtained for 
gamma-ray detection (typically 200 /is) and time tagging 
with uncertainty near ^1 fis. The Super AGILE hard X- 
ray monitor is positioned on top of the GRID. Super- 
AGILE is a coded aperture instrument operating in the 
18-60 keV energy band with about 15 mCrab sensitivity 
in one day integration, 6 arcmin angular resolution and 
~1 sr field of view (Feroci et al. 2007; Costa et al. 2001). 

In this work we analyze all available AGILE /GKID 
data suitable for timing analysis collected up to 2008 
April 10 for the four known gamma-ray pulsars included 
in the AGILE Team source list:^^ Vela, Crab, Geminga 
and B 1706 44. The other two EGRET pulsars, B 1055 
52 and B 1951+32, are part of the AGILE Guest Ob- 
server program. As expected, only the Crab pulsar has 
been detected by SuperAGILE and the X-ray data have 
been used to cross-check and test AGILE timing perfor- 
mances. 

The AGILE observations are presented in § 2, as well 
as the criteria for photons selections. The observations 
and the timing analysis from the parallel radio and X-ray 

See http://agile.asdc.asi.it for details about AGILE Data Pol- 
icy and Target List. 



observations of the four targets are described in § 3. The 
procedures for the timing analysis of gamma-ray data 
are introduced in § 4 and the results of their application 
are reported in § 5, where timing calibration tests are 
also dealt with. Discussion of the scientific results and 
conclusions are the subjects of §§ 6 and 7, respectively. 

2. AGILE OBSERVATIONS AND DATA 
REDUCTION 

The AGILE spacecraft was placed in a Low Earth Or- 
bit (LEO) at ^535 km mean altitude with inclination 
^2.5°. Therefore, Earth occultation strongly affects ex- 
posure along the orbital plane, as well as high particle 
background rate during South Atlantic Anomaly tran- 
sits. However, the exposiue efficiency is >50% for most 
AGILE revolutions. AGILE pointings consist of long 
exposures (typically lasting 10 30 days) slightly drifting 
(<l°/day) with respect to the starting pointing direction 
in order to match solar-panels illumination constraints. 
The relatively uniform values of the effective area and 
point spread function within ~40° from the center of the 
field of view of the GRID, allow for one- month point- 
ings without significant vignetting in the exposure of the 
target region. 

Pulsar data were collected during the mission Science 
Verification Phase (SVP, 2007 July-November) and early 
pointings^^ (2007 December 2008 April) of the AO 1 Ob- 
serving Program. It is worth noting that a single AGILE 
pointing on the Galactic Plane embraces about one-third 
of it, allowing for simultaneous multiple source targeting 
(e.g. the Vela and Anti-Center regions in the same field 
of view with Crab, Geminga and Vela being observed at 
once; Figure 1). 

The AGILE Commissioning and Science Verification 
Phases lasted about seven months from 2007 April 23 
to November 30 including also Instrument Time Calibra- 
tion. On 2007 December 1 baseline nominal observations 
and a pointing plan started together with the Guest Ob- 
server program AOl. Timing observations suitable for 
pulsed signal analysis of the Vela pulsar started in mid 
2007 July (at orbit 1,146) after engineering tests on the 
pay load. 

The Vela region was observed (with optimal exposure 
efficiency) for ^40 days during the SVP and again for 
~30 days in AO 1 pointing number 3 (2008 January 8- 
February 1). PSR B 1706-44 was within the Vela and 
Galactic Center pointings for ^30 days during the SVP 
and for ^45 days during AO 1 pointings numbers 5 and 
6 (2008 February 14-March 30). The Anti-center re- 
gion (including Crab and Geminga) was observed for ~40 
days, mostly in 2007 September, and in 2008 April (AO 1 
pointing number 7) with the addition of other sparse 
short Crab pointings for SuperAGILE calibration pur- 
poses during the SVP (see Table 1 and Figure 2 for de- 
tails about targets coverage). 

GRID data of the relevant observing periods were 
grouped in 20 subsets of 200 orbits each (correspond- 
ing to ~15 days of observation) starting from orbit 
1,146 (54,294 MJD). Data screening, particle background 

2'^ 1. Cygnus Field 1 (Z = 89°, = 9.9°), 2. Virgo Field (l = 
264°, f) = 56.5°), 3. Vela Field (« = 283°, fe = -6.8°), 4. South Gal. 
Pole (; = 240°, 6 = -50°), 5. Musca Field {I = 303°, b = -9°), 
6. Gal. Center (( = 332°, b = 0°), 7. Anti-Center (I = 193°, 
6 = 8.1°). 
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Fig. 1. — Gaussian-smoothed AGILE intensity map (~120° X 60°, units: ph cm s sr , E>100 MeV) in Galactic coordinates 
integrated over the whole observing period (2007 July 13-2008 April 10) and centerd at I = 223°, b = 0°. The AGILE field of view (radius 
~60°) can embrace in a single pointing Vela (l = 263.6°, b = -2.8°), Geminga (I = 195.1°, b = 4.3°) and Crab (/ = 184.6°, b = -5.8°) as 
well as diffuse emission from the Galactic Disk. 



PRE-WHITENING 



E 




5420D 54300 54400 54500 
Date (MJD) 



54200 54300 54400 54500 
Date (MJD) 



Fig. 2. — Post-fit timing residuals (in milli-turns) as a function of 
the Modified Julian Day resulting from the observation at 1.4 GHz 
of the three radio pulsars which are discussed in this paper: from 
the top Vela, Crab and PSR B1706-44. The panels in the left 
(right) column report the residuals of the best available timing 
solutions obtained over the data span not including (including) 
the correction of the timing noise via the use of the AR term (see 
§ 4). Note that for the Crab pulsar the scale on the vertical axis 
of the panel in the right column is amplified by a factor 10 with 
respect to that in the left panel. The time intervals corresponding 
to the useful AGILE pointings for each target are also given as the 
black sections of the bar at the bottom of each panel. 



filtering and event direction and energy reconstruc- 
tion were performed by the AGILE Standard Analysis 
Pipeline (BUILD-15) for each subset with an exposure 
>10^ cm^ s at i? > 100 MeV. Observations affected by 
coarse pointing, non-nominal settings or intense parti- 
cle background (e.g. orbital passages in South Atlantic 
Anomaly) and albedo events from the Earth's limb were 
excluded from the processing. 

A specific optimization on the events extraction pa- 
rameters is performed for each target in order to max- 
imize the signal to noise ratio for a pulsed signal. The 
optimal event extraction radius around pulsar positions 
varies as a function of photon energy (and then it is re- 
lated to pulsar spectra) according to the Point Spread 
Function. However, for E > 100 MeV broad band anal- 
ysis, a fixed extraction radius of ~5° (a value slightly 
higher than the Point Spread Function 68% contain- 
ment radius) produces comparable results with respect 
to energy-dependent extraction. 

Quality flags define different GRID event classes. The 
G event class includes events identified with good confi- 
dence as photons. Such selection criteria correspond to 
an effective area of ^250 cm^ above 100 MeV (for sources 
within 30 deg from the center of the field-of-view). The 
L event class includes events typically affected by an or- 
der of magnitude higher particle contamination than G, 
but yielding an effective area of ~500 cm^ at E > 100 
MeV, if grouped with the G class. We performed our 
timing analysis looking for pulsed signals using both G 
class events and the combination of G-l-L events. In gen- 
eral, the signal-to-noise ratio of the pulsed signal is maxi- 
mized using photons collected within ^^40° from the cen- 
ter of GRID field of view and selecting the event class 
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TABLE 1 

Observation parameters of the relevant data subsets (grouped in uniterrupted 

TARGET observations) FOR EACH ANALYZED PULSAR CONSIDERING FOR COMPARISON ALL 

G+L CLASS EVENTS (SEE TEXT) WITH E > 100 MeV OBSERVED WITHIN 60° FROM THE 
CENTER OF THE FIELD OF VIEW AND EXTRACTED WITHIN 5° FROM THE PULSAR POSITION. 



PSR 


ObsID* 


-L FI RST 


-L LAST 


\^/ 


Total Counts'^ 


Exposure^ 






[iVl J UJ 


ri\,f T"ni 

[iVi J UJ 


[deg] 




[lu cm sj 


Vela 


SVP 1 


54,294.552 


54,305.510 


26.2 


6440 


1.87 


Vela 


SVP 2-3-4 


54,311.512 


54,344,569 


40.9 


12651 


3.26 


Vela 


SVP 5 


54,358.525 


54,359.521 


41.4 


391 


0.11 


Vela 


SVP 6 


54,367.525 


54,377.521 


48.9 


1355 


0.26 


Vela 


SVP 8 


54,395.554 


54,406.501 


55.3 


113 


0.03 


Vela 


AOl 2 


54,450.540 


54,457.157 


56.2 


564 


0.01 


Vela 


AO 1 2-3 


54,464.784 


54,480.966 


22.1 


4726 


3.13 


Vela 


AO 1 3-4 


54,492.948 


54,505.377 


43.2 


6293 


1.16 


Vela 


AO 1 4-5-6 


54,508.528 


54,528.537 


43.6 


9747 


2.22 


Vela 


AOl 6 


54,546.095 


54,.561.427 


41.7 


703 


0.17 


Vela 


total° 


54,294.552 


54,-561. 127 


36.1 


42983 


12.22 


Vela 


gamma' 


54,294.552 


54,561.427 


36.1 


6140 


5.28 


Geminga 


SVP 2 


54,308.871 


54,314.507 


55.0 


187 


0.06 


Geminga 


SVP 3 


54,324.514 


54,335.508 


41.3 


781 


0.24 


Geminga 


SVP 4^5-6-7 


54,344.518 


54,351.235 


22.2 


17823 


6.19 


Geminga 


SVP 8 


54,395.520 


54,406.504 


37.6 


1926 


0.48 


Geminga 


AOl 6 


54,528.531 


54,529.270 


50.7 


38 


0.01 


Geminga 


AOl 6 


54,546.093 


54,549.428 


39.0 


132 


0.03 


Geminga 


AO 1 6-7 


54,555.514 


54,566.5 


11.7 


4860 


2.03 


Geminga 


totaP 


54,308.871 


54,566.5 


20.2 


25962 


9.04 


Geminga 


gamma*^ 


54,308.871 


54,566.5 


20.2 


3874 


4.54 


Crab 


SVP 1-2 


54,305.509 


54,314.50 


50.2 


2404 


0.58 


Crab 


SVP 3 


54,324.512 


54,335.508 


28.2 


761 


.0.32 


Crab 


SVP 4-5-6-7 


54,344.518 


54,386.502 


23.9 


20499 


6.41 


Crab 


SVP 8 


54,395.518 


54,406.506 


45.4 


1537 


0.32 


Crab 


AO 1 3-4 


54,494.447 


54,505.390 


55.8 


350 


0.05 


Crab 


AOl 4 


54,508.506 


54,510.474 


54.7 


359 


0.06 


Crab 


AOl 6 


54,528.531 


54,529.140 


46.9 


100 


0.02 


Crab 


AOl 6 


54,546.091 


54,549.423 


41.5 


176 


0.03 


Crab 


AO 1 6-7 


54,555.516 


54,566.5 


21.8 


4973 


1.89 


Crab 


total° 


54,305.509 


54,566.5 


26.1 


31159 


9.68 


Crab 


gamma' 


54,305.509 


54,566.5 


26.1 


4062 


4.11 


B 1706-44 


SVP 1 


54,294.552 


54,305.503 


50.9 


4641 


0.91 


B1706-44 


SVP 2-3-4 


54,311.531 


54,322.891 


34.5 


15657 


4.39 


B1706-44 


SVP 7-8 


54,386.524 


54,393.752 


41.5 


5521 


1.28 


B1706-44 


AO 1 2-3 


54,470.292 


54,480.965 


53.7 


2169 


1.46 


B1706-44 


AOl 3 


54,492.948 


54,497.509 


44.9 


1507 


0.20 


B1706-44 


AO 1 5-6 


54,510.482 


54,555.509 


22.4 


24219 


7.36 


B1706-44 


total" 


54,294.552 


54,555.509 


28.7 


53714 


15.6 


B1706-44 


gamma' 


54,294.552 


54,555.509 


28.7 


8463 


6.56 



^Observation ID: SVP=Science Verification Phase (grouped in subset of 200 orbit each 
starting from 1,146), AO l=Scientific Observations Program pointings (see AGILE Mission 
Announcement of Opportunity Cycle-1; http;//agile. asdc.asi.it/). 

''Mean off-axis angle. 

'^Source photons + diffuse emission photons + particle background. 
''Good observing time after dead-time and occultatioii corrections. 

"^Total G+L class events with E > 100 MeV, 5° max from pulsar position, 60° max from 
field of view center. 

'High-confidence photon events (G class only) with E > 100 MeV, 5° max from pulsar 
position, 60° max from field of view center. 



G. For very strong sources, such as Vela, or sources lo- 
cated in low background regions, it is possible to include 
also photons in the 40-60° off-axis range and belonging 
to the G-hL event class typically improving the count 
statistics by up to a factor of three (obviously imply- 
ing also a much higher background), without affecting 
or even improving the detection significance. For each 
pulsar, Table 1 summarizes the exposure parameters of 
the relevant data subsets grouped in uninterrupted tar- 
get observations. For simplicity and comparison, we re- 
ported in Table 1 and Table 2 observation and detection 
parameters obtained with the same extraction criteria 
and energy range. For all observations with the targets 



within 60° from the pointing direction, we list the av- 
erage angular distance (9) of the source position from 
the pointing direction, the sum of all G-l-L events with 
E > 100 MeV whose direction is within 5° from the 
pulsar position and the G only events ("gamma pho- 
tons"), selected with the same criteria. The choice to 
include all G-l-L events yields a "dirty" data set: for the 
case of Vela, ^10% of the total counts are ascribable to 
the Galactic Plane, and >50% are particle background 
(~10,000 source counts). The contamination is reduced 
in the "gamma" entry, characterized by ^20% diffuse 
emission and ^5% particle background (~5,000 source 
counts). The corresponding exposure for each data sub- 
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set was calculated with the GRID scientific analysis task 
AG_ExpmapGcn according with the above parameters. 
Summing up the photon numbers, we see that the over- 
all photon statistics accumulated (accounting for back- 
ground contamination) is comparable to that collected 
by EGRET for the same four pulsars. 

A maximum likelihood analysis (ALIKE task) on the 
AGILE data from the sky areas containing the four 
gamma-ray pulsars yielded source positions, fluxes and 
spectra in good agreement with those reported both in 
the EGRET catalogue for the corresponding 3EG sources 
(Hartman et al. 1999) and in the revised catalogue by 
Casandjian & Grenier (2008). As an example, Fig- 
ure 1 reports an AGILE intensity map displaying Vela, 
Geminga and Crab. Details on gamma-ray imaging and 
spectra of the four sources will be reported in future pa- 
pers when count statistics significantly higher than that 
currently available are collected and in-flight Calibration 
files are finalized. 

In this paper we focus on timing analysis. However, it 
is worth noticing that pulsed counts provide a gamma- 
ray pulsar flux estimate independent from the likelihood 
analysis, as described in § 5. 

3. RADIO/X-RAYS OBSERVATIONS AND TIMING 

In order to perform A GILE timing calibration through 

accurate folding and phasing, as described in §§ 4 and 5, 
pulsar timing solutions valid for the epoch of AGILE ob- 
servations were required. Thus, a dedicated pulsar mon- 
itoring campaign (that will continue during the whole 
AGILE mission) was undertaken, using two telescopes 
(namely Jodrell Bank and Nancay) of the European Pul- 
sar Timing Array (EPTA), as well as the Parkes radio 
telescope of the Commonwealth Scientific and Industrial 
Research Organisation (CSIRO) Australia Telescope Na- 
tional Facility (ATNF) and the 26m Mt Pleasant radio 
telescope operated by the University of Tasmania. 

In particular, the observations of the Vela Pulsar have 
been secured by the Mt. Pleasant Radio Observatory. 
Data were collected at a central frequency of 1.4 GHz. 
Given the pulsar brightness, it is possible to extract a 
pulse Time of Arrival (ToA) about every 10 seconds, so 
that a total of 4,098 ToAs have been obtained for the 
time interval between 2007 February 26 (MJD 54,157) 
and 2008 March 23 (MJD 54,548), encompassing the 
whole time span of the AGILE observations. During this 
time interval the pulsar experienced a small glitch (frac- 
tional frequency increment Aug/u ~ 1.3 x 10^^), which 
presumably happened around 2007 August 1 (with an 
uncertainty of ±3 days due to the lack of observations 
around this period) (see § 6.3). Ephemeris for the pre- 
and post-glitch time intervals were then separately cal- 
culated and applied to the gamma-ray folding procedure. 

The observations of the Crab Pulsar have been pro- 
vided by the telescopes of Jodrell Bank and Nangay. At 
Jodrell Bank the Crab Pulsar has been observed daily 
from 2006 December 09 (MJD 54,078) to 2007 Octo- 
ber 06 (MJD 54,379) and again from 2008 February 06 
(MJD 54,502) to 2008 April 10 (MJD 54,566) which re- 
sult in 334 ToAs. The observations were mainly per- 
formed at 1.4 GHz with the 76 m Lovell telescope, with 
some data also taken with the 12-m telescope at a cen- 
tral frequency of 600 MHz. The Nangay Radio Telescope 
(NRT) is a transit telescope with the equivalent area of a 



93 m dish and observed at a central frequency of 1.4 GHz 
over the time interval between 2006 December 05 (MJD 
54,074) and 2007 September 07 (MJD 54,350), produc- 
ing 64 ToAs. A timing solution was produced by joining 
the ToAs from the two telescopes and accounting for the 
phase shift that naturally ensues from data sets coming 
from different telescopes. 

The observations of pulsar B 1706—44 have been per- 
formed at the 64 m telescope at Parkes, in Australia, at 
the mean observing frequency of 1.4 GHz. They pro- 
duced a total of 20 ToAs over the interval from 2007 
April 30 (MJD 54,220) to 2008 Apr 6 (MJD 54,562), 
thus covering the whole A GILE observing time span for 
this target. 

The timing of all pulsars is performed using the 
TEMP02 software (Hobbs et al. 2006; Edwards et al. 
2006). It first converts the topocentric TOAs to solar- 
system barycentric TOAs at infinite frequency^'^ (using 
the Jet Propulsion Laboratory [JPL] DE405 solar-system 
ephemeris^^) and then performs a multi-parameter least- 
square fit to determine the pulsar parameters. The differ- 
ences between the observed barycentric ToAs and those 
estimated by the adopted timing model are represented 
by the so-called residuals. The procedure is iterative 
and improves with longer timescales of observations. An 
important feature developed by TEMP02 is the possi- 
bilty of accounting for the timing noise in the fitting 
procedure (see § 4). This is particularly useful in tim- 
ing the young pulsars: in fact most of them suffer of 
quasi-random fluctuations (typically characterized by a 
very red noise spectrum) in the rotational parameters, 
whose origin is still debated. TEMP02 corrects the ef- 
fects of the timing noise on the residuals by modeling its 
behaviour as a sum of harmonically correlated sinusoidal 
waves (see Equation [2] in § 4) that is subtracted from 
the residuals. 

The left panels (labeled as pre-whitening) in Figure 2 
reports the timing residuals obtained over the data span 
of the radio observations without the correction for 
the timing noise, whereas the right panels (labeled as 
whitened) show the residuals after the application of the 
correction. The comparison between the panels in the 
two columns shows that this procedure has been very ef- 
fective in removing the timing noise in Vela, Crab and 
PSR B1706— 44. The impact of this whitening procedure 
on the gamma-ray data analysis is discussed in § 5. 

Geminga is a radio-quiet pulsar whose ephemeris can 
be obtained from X-ray data. Following the demise of 
CGRO, Geminga was regularly observed with XMM- 
Newton in order to maintain its ephemeris for use in 
analyzing observations at other wavelengths. We ana- 
lyzed all eight observations of Geminga (IE 0630+178, 
PSR J 0633+1746) taken with the X-ray (0.1 15 keV) 
cameras of XMM-Newton (Jansen et al. 2001) between 
MJD 52,368 and MJD 54,534 (2002 April - 2008 March), 
with exposures in the range 20-100 ks. 

The data were processed using version 7.1.0 of the 
XMM-Newton Science Analysis Software (SAS) and the 
calibration files released in 2007 August. For the timing 

Dispersion measure is obtained as part of timing solution 
for Crab (DM=56.76(1)), Vela (DM=68.15(2)) and from Johnston 
et al. (1995) for B 1706-44 (DM=75.69(5)). 

'^'^ See ftp;//ssd.jpl. nasa.gov/pub/eph/export/DE405/de405.iom/. 
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TABLE 2 

The photon harvest from the AGILE observations of gamma-ray pulsars. 



PSR P E d Pulsed Counts'' Xr (d.o.f) Exposure'' Pulsed Flux<= 

(ms) (erg s"!) (kpc) (10* cm^ s) 10"* ph. cm-^ s"! 

Vela 89.3 6.92 x lO^^ 0.29 9, 170 ± 580 225.51 (9) 12.22 940 ± 60 

Geminga 237.1 3.25 x 103" 0.16 2, 200 ±480 10.44 (9) 9.04 300 ± 70 

Crab 33.1 4.61x10'''* 2.00 2, 120 ± 5.30 10.71 (9) 9.68 270 ± 70 

B1706-44 102.5 3.41 x 10^" 1.82 2, 370 ± 720 9.11 (9) 15.6 190 ± 60 

Note. — See text for details on data reduction and timing analysis. 

^'Pulsed counts (G+L event class) with E > 100 MeV, 5° max from pulsar position, 60° max from field of 
view center, 10 bins. 
''Good observing time after dead-time and occultation corrections. 

"^Calculated with the expression Cpf/E, where Cp=pulsed counts, E=exposure, /=factor accounting for 
source counts at angular distance >5° from source position according to the point spread function (/ ~ 1.25). 



analysis wc could only use the pn data (operated in Small 
Window mode: time resolution 6 ms, imaging on a 4' x 4' 
field), owing to the inadequate time resolution of tlie 
MOS data. We selected only single and double photon 
events (patterns 0-4) and applied standard data screen- 
ing criteria. Photons arrival times were converted to the 
Solar system barycenter using the SAS task barycen. To 
extract the source photons we selected a circle of 30" ra- 
dius, containing about 85% of the source counts. Using 
standard folding and phase- fitting techniques, the source 
pulsations were clearly detected in all the observations. 
We derived for Geminga the long-term spin parameters 
valid for the epoch range 52,369-54,534 MJD. Note that 
the absolute accuracy of the XMM-Newton clock is bet- 
ter than 600 ^s (Kirsch et al. 2004). 

The Crab Pulsar is embedded in the Crab Nebula and 
represents about 10% of its flux in the hard X-ray band. 
Being a bright source with a relatively soft high-energy 
spectrum, the Crab pulsar is easily detected by Super- 
AGILE in less than one AGILE orbital revolution. Since 
the on-board time reconstruction of SuperAGILE (Feroci 
et al. 2007) is different from that of the GRID, as a cross- 
check and test of the SuperAGILE timing performances 
we processed the X-ray data of an on-axis observation 
of ~0.7 days (54,360.7-54,361.4 MJD), corresponding to 
an effective exposure of ^^41 ks. In the analysed data 
the passages of AGILE in the South Atlantic Anomaly 
and intervals of source occultation by the Earth were 
excluded. The time entries in SuperAGILE event files 
are in the on-board time reference system (Coordinated 
Universal Time) and were converted in Terrestrial Dy- 
namical Time before being processed and analysed by 
the same procedures used for the GRID (see § 4). A 
complete analysis of the X-ray observations of Crab with 
SuperAGILE is beyond the scope of this work: it will be 
presented in a future paper, as well as the search for the 
Crab pulsed signal in the AGILE mini-calorimeter data. 

4. GAMMA-RAY TIMING PROCEDURES 

In this section we will describe the timing procedures 

we have adopted. They have been implemented with two 
aims: to verify the timing performances of AGILE mid to 
maximize the quality of the detection of the four targets. 

A GILE on-board time is synchronized to Coordinated 
Universal Time (UTC) by Global Positioning System 
(GPS) time sampled at a rate of 1 Hz. Arrival time 
entries in AGILE event lists files are then corrected to 
Terrestrial Dynamical Time (TDT) reference system at 



ground segment level. In order to perform timing analy- 
sis, they have also to be converted to Barycentric Dynam- 
ical Time (TDB) reference and corrected for arrival de- 
lays at Solar System Barycenter (SSB). This conversion is 
based on the precise knowledge of the spacecraft position 
in the Solar System frame. To match instrumental mi- 
crosecond absolute timing resolution level, the required 
spacecraft positioning precision is <0.3 km. This goal is 
achieved by the interpolation of GPS position samples 
extracted from telemetry packets. Earth position and 
velocity with respect to SSB are then calculated by JPL 
planetary ephemeris DE405. All the above barycentric 
corrections are handled by a dedicated program (imple- 
mented in the AGILE standard data reduction pipeline) 
on the event list extracted according to the criteria de- 
scribed in § 2. 

Pulsed signals in GRID gamma-ray data cannot be 
simply found by Fourier analysis of the photon SSB ar- 
rival times, since the pulsar rotation frequencies are 4-5 
order of magnitude higher than the gamma-ray pulsars 
typical coimt-rates (10 100 coimts/day). The determi- 
nation of the pulsar rotational parameters in gamma-ray 
must then start from an at least approximate knowledge 
of the pulsar spin ephemeris, provided by observations at 
other wavelengths. 

Standard epoch folding is performed over a tri- 
dimensional grid centered on the nomimal values of the 
pulsar spin frequency and of its first and second or- 
der time derivatives, uq and z>b, as given by the as- 
sumed (radio or X-ray) ephemeris at their reference 
epoch to = PEPOCH. The axes of the grid are explored 
with steps equal to l/Tgpan, '^/T^%an ^/"^Ipan respec- 
tively (all of them oversampled by a factor 20), where 
T^span is the time span of the gamma-ray data. For any 
assigned tern [v; v\ i>], the pulsar phase associated to 
each gamma-ray photon is determined by the expression: 

$* = $0 -I- vM -h -!>At^ -i>Af^ (1) 
2 6 

where At = t — to is the difference between the SSB 

arrival time t of the photon and the reference epoch 
to of the ephemeris and $o is a reference phase (held 
fixed for all the set). A light-curve is formed by bin- 
ning the pulsar phases of all the photons and plotting 
them in a histogram. Pearson's statistic is then ap- 
plied to the light-curves resulting from each set of spin 
parameters, yielding the probabilities of sampling a uni- 
form distribution. These probabilities V{y\ v; i)) are 
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then weighted for the number AQt of steps over the grid 
which has been necessary to reach the set [v; i)] i)] start- 
ing from [i/q; uq; Vq]. The maximum value over the grid 
of iS = 1 — MstP{y; v\ i>) finally determines which are 
the best gamma-ray rotational parameters for the target 
source in the surroundings of the given ephemeris. Of 
course, the higher the value of 5, the higher is the statis- 
tical significance of a pulsating signal. We note that this 
approach allows to avoid any arbitrariness in the choice 
of the range of the parameters to be explored, which oth- 
erwise can affect the significance of a detection. 

Due to the brightness of the sources discussed in this 
paper, period folding around the extrapolated pulsar spin 
parameters obtained from publicly available ephemerides 
(ATNF Pulsar Catalog,3° Manchester et al. 2005, and Jo- 
drell Bank Crab ephemeris archive, Lyne et al. 1993) 
or from recent literature led us to firmly detect (>5(t) 
the pulsations for all the four targets with a reasonable 
number of trials (^1,000). However, for three of them 
(Crab, Vela and PSR B1706-44), the best gamma-ray 
spin period fell outside the 2>a uncertainty range of the 
adopted radio ephemeris, and the detection significance 
S resulted lower than that derived using contemporary 
timing solutions (accounting for the elfects of timing 
noise and/or occasional glitches), provided by dedicated 
radio observation campaigns (see § 3). As expected, for 
steady and older pulsars, such as Geminga, the availabil- 
ity of contemporary rotational parameters turned out to 
be less important; even using a few-year-old ephemeris, 
the pulsed signal could be detected within only <100 pe- 
riod search trials around the extrapolated X-ray timing 
solutions. 

An additional significant improvement (see § 5) in the 
detection significance has been obtained by accounting 
also for the pulsar timing noise in the folding procedure. 
This exploits a tool of TEMP02 (Hobbs et al. 2006; Ed- 
wards et al. 2006), namely the possibility of fitting tim- 
ing residuals with a polynomial harmonic function Ai? 
in addition to standard positional, rotational and (when 
appropriate) binary parameters (Hobbs et al. 2004): 

N 

AR{At) = ^ flfe sin(A:a;At) + bk cos{kujAt) (2) 

where A'^ is the number of harmonics (constrained by 
precision requirements on radio timing residuals, as well 
as by the span and the rate of the radio observations), 
Gk and bk are the fit parameters (i.e. the WAVE terms 
in TEMP02 ephemeris files), and u> = 27r(Tradio(l + 
4/A^))~-'^ is the main frequency (i.e. WAVEqm in 
TEMP02 ephemeris files) related to the radio data time- 
span Tradio- 

If the spin behavior of the target is suitably sampled, 
the harmonic function AR can absorb the rotational ir- 
regularities of the source, in a range of timescales ranging 
from ^Tradio down to about the typical interval between 
radio observations. As an example, the peak-to-peak am- 
plitude of the AR fluctuations for Crab related to the 
radio monitoring epochs (54,074-54,563 MJD) covering 
our AGILE observations are of the order of ^1 ms, corre- 
sponding to a phase smearing >0.03, a value significantly 

See http;//www.atnf.csiro.au/research/pulsar/psrcat/. 
See http://www.jb.man.ac.Uk/~pulsaj:/crab.html. 



affecting the time resolution of a >50 bins light-curve. 
Under the assumption that the times of arrival of the 
gamma-ray photons are affected by timing noise like in 
the radio band, gamma-rays folding can properly account 
for AR, extending the relation 1 to: 

^ = ^* + {v + vAt + h)At^)AR 

~ $0 + i^At{l + ^) + \vAe + \i>At^ (3) 
At I 

As reported in § 5, this innovative phasing technique 
improves significantly the gamma-ray folding accuracy 
for young and energetic pulsars, especially when using 
long data spans, like those of the AGILE observations. 
Of course, the implementation of this procedure requires 
radio observations covering the time span of the gamma- 
ray observations making the radio monitoring described 
in § 3 all the more important. 

5. TIMING CALIBRATION TESTS AND RESULTS 

In order to verify the performances of the timing anal- 
ysis procedure described in § 4, a crucial parameter to 
check is the difference between pulsar rotation param- 
eters derived from radio. X-ray and gamma-ray data. 
Figure 3 shows the AGILE period search result for the 
Crab pulsar (corresponding to the MJD 54,305-54,406 
observations, significantly affected by timing noise, as 
shown in Figure 2). The implementation of the fold- 
ing method described in section § 4 (including timing 
noise corrections as given by Equation [3]) allowed for 
a perfect match between the best period resulting from 
gamma-ray data and the period predicted by the radio 
ephemeris with discrepancies APcrab ~ 3 x 10^^^ s, com- 
parable to the period search resolution rcrab ~ 2 x 10"-^^ 
s. This represents also an ultimate test for the accu- 
racy of the on-board AGILE Processing and Data Han- 
dling Unit (PDHU; Argan et al. 2004) time management 
(clock stability in particular) and on-ground barycen- 
tric time correction procedure. Standard folding with- 
out AR term implies radio-gamma period discrepancies 
one order of magnitude higher (APcrab ~ 4 x 10~^^ 
s). Moreover, it lowers also the statistical significance of 
the detection and effective time resolution of the light- 
curve. For the Crab, the value of the Pearson's 
statistic introduced in § 4 (we quote reduced val- 
ues) goes from ~6.3 (when using the AR term) down 
to ^^4.2 (ignoring the AR term) when folding the data 
into a 50 bins light-curve (see Figure 4). Obviously, ig- 
noring timing noise in the folding process would yield 
discrepancies (and light curve smearing) which are ex- 
pected to grow when considering longer observing time 
span. Thus the contribution of timing noise should be 
considered both in high-resolution timing analysis and in 
searching for new gamma-ray pulsars. The same analy- 
sis applied to Vela and PSR B 1706-44 - much less af- 
fected than Crab by timing noise in the considered data 
span - led to similar results for the period discrepancies 
(APvcia 8 X 10^12 APbi706 ~ 3 X 10"", to be com- 
pared with the period search resolution rveia ~ 9 x 10~^^ 
s and rBi706 ~ 2 x 10~^^ s). 

For the radio-quiet Geminga pulsar we used X-ray 
ephemeris obtained from XMM-Newton data (see § 3) as 
a starting point for the period search. Due to the stabil- 
ity of the spin parameters of this relatively old pulsar, not 
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P_trial-P_radio [s] 

Fig. 3. — Period search result for the Crab pulsar (period trials 
vs. Pearson statistics). The radio period (vertical line at Ptrial~ 
^Vadio = 0) is 33.607554009(4) ms {PEPOCH = 54,362.242 
MJD). Dashed line is obtained from standard folding period tri- 
als, while continuos line is from the folding technique accounting for 
timing noise (see § 3). The new method allow us to perfectly match 



the radio period in gamma rays (P^ 



gamma_tn 



3 X 10 



-12 



s). It is worth noting that in the considered data span (54,305- 
54,406 MJD) the period search resolution (rc^ab ~ 2 X 10^^"^ s) is 
about an order of magnitude higher than the Icr error on the radio 
period. 




Fig. 4. — The Crab light-curve (50 bin) corresponding to the 
data span 54,305-54,406 MJD obtained by folding including AR 
terms compared with that obtained neglecting timing noise (dashed 
line). In observing periods strongly ailected by timing noise, the 
smearing effects reduce detection significance and observed pulsed 
counts (see text). 



significantly affected by timing noise, WAVE parameters 

are not required for the folding process. The X-ray versus 
gamma-ray period discrepancy was APccm ^ 9 x 10~^^ 
s, whereas rQcm 7x 10~^^ s. Wc hero note that the fre- 
quency resolution l/Tgpan ~ 10~^ Hz of AGILE is about 
one order of magnitude better than that corresponding 
to a single XMM-Newton exposure (lasting < 100 ks), 
but the few year long X-ray data span implies an overall 
much better effective resolution 1/Txmm ~ 5 x 10~^ Hz 
for the XMM-Newton data. 

The resulting gamma-ray light-curves, covering differ- 
ent energy ranges for the four pulsars, are shown in Fig- 
ures 5, 6, 7, and 8. The pulsed flux was computed con- 
sidering all the counts above the minimum of the light- 




FlG. 5. — Vela pulsar light-curves (P ~ 89.3 ms) for different en- 
ergy bands {E < 100 McV, 40 bins, resolution: ~2.2 ms; E > 100 
MeV, 100 bins, resolution: ~0.9 ms; E > 1 GcV, 20 bins, resolu- 
tion: ^4.5 ms; G+L class events) obtained by integrating all avail- 
able post-glitch data (54,320-54,561 MJD). The radio cphcmcris 
and the 8,192 bins light-curve (bottom panel) arc obtained by the 
analysis of ~4,100 ToAs observed at Mt. Pleasant Radio Observa- 
tory in Tasmania (radio observation interval 54,157-54,548 MJD). 



curve, using the expression PF = Ctot 



-nN„ 



sociated error app = (Ctot +" ^Ar„i„) — "(-^min) 



niin and its as- 

0^ n(iV„,i„)i/2, 

where Ctot arc the total counts, n is the number of bins 
in the light-curve and A^min are the counts of bin cor- 
responding to the minimum. This method is "bin de- 
pendent", but reasonable different choices of both the 
number of bins (i.e. n > 10) and of the location of the 
bin center (10 trial values were explored for each choice 
of n) do not significantly affect the results. Several mod- 
els (polar cap, slot gap) predict that gamma-ray pulsar 
emission is present at all phases. When enough counts 
statistics will be available, (e.g. ~10,000 counts for the 
Crab pulsar), it will be possible to estimate the unpulsed 
emission (due to the pulsar plus a possible contribution 
from the pulsar wind nebula) by considering the differ- 
ence among the total source flux obtained by a likelihood 
analysis on the images and the pulsed flux estimated with 
the above method. 

Pulsed counts and related Pearson statistics for the 
four pulsars are reported in Table 2 (for the standard 
event extraction parameters as in Table 1). The result- 
ing fluxes (Pulsed Counts/Exposure) are consistent with 
those reported in the EGRET Catalogue (Hartman et al. 
1999). We note that source-specific extraction param- 
eters (event class, source position in the field of view. 
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Fig. 6. — Geminga light-curves (P ~ 237.1 ms) for different 
energy bands {E < 100 McV, 20 bins, resolution: ^11. 8 ms; E > 
100 MeV, 100 bins, resolution: ~2.4 ms; B > 1 GcV, 20 bins, 
resolution: ~11.8 ms; G class events). The X-ray eplicmcris and 
the 1—8 keV 40 bins light-curve (bottom panel) are obtained by the 
analysis of XMM-Newton data (observation interval 52,369-54,534 
MJD). 



energy band etc.) which maximize reduced can sig- 
nificantly improve detection significance. For example, 
including only event class G in the timing analysis of 
Crab and Geminga halves the number of pulsed counts 
while doubling the values. 

Despite the verj^ satisfactory matching of the pulsar 
spin parameters found in radio (or X-ray) and gamma- 
ray (supporting the clock stability and the correctness 
of the SSB transformations), possible systematic time 
shifts in the AGILE event lists could be in principle af- 
fecting phasing and must be checked. For example, an 
hypothetical constant discrepancy of terr of the on-board 
time with respect to UTC would result in a phase shift 
^err = [terr mod P)/P, where P is the pulsar period. 
The availability of radio observations bracketing the time 
span of the gamma-ray observations (or of X-ray ob- 
servations very close to the gamma-ray observations for 
the case of Geminga) allowed us to also perform accu- 
rate phasing of multi-wavelength light-curves. In doing 
that, radio ephemeris reference epochs were set to the 
main peak of radio light-curves at phase ^peak = 0. In 
view of Equations (2) and (3), this is achieved by setting 

$0 = -i/X;f=i H (typically $o < 10"^). We found that 

the phasing of the AGILE light-curves of the four pul- 
sars (radio/X-r ays/gamma-ray peaks phase separations) 
is consistent with EGRET measurements (Fierro et al. 



Fig. 7. — Crab pulsar light-curves (P ~ 33.1 ms) for different 
energy bands {E > 100 McV, 50 bins, resolution: ~0.7 ms; E > 500 
McV, 40 bins, resolution: ~0.8 ms; G class events) obtained by 
integrating all available data presented in Tabic 1. The X-ray (I8- 
60 keV) SuperAGILE 50 bins light-curve is obtained from a ~41 ks 
observation taken at ~54,361 MJD. The radio eplicmcris and the 
2048 bins light-curve (bottom panel) are obtained by the analysis 
of 334 ToAs observed by Jodrell Bank and Nangay radio telescopes 
(observation interval 54,078-54,566 MJD). 

1998; Thompson et al. 1996; Jackson & Halpern 2005, 
see § 6 for details) implying no evidence of systematic 
errors in absolute timing with an upper limit terr < 1 
ms. 

Comparison with the SuperAGILE light-curve peak 
(see § 3) is also interesting (SuperAGILE on-board time 
processing is more complex than that of the GRID). The 
Crab SuperAGILE light-curve (Figure 7) was produced 
with the same folding method reported in § 4, yielding 
63, 700±8, 700 pulsed counts (~3% of the total counts in- 
cluding background). Inspection of Figure 7 shows that 
the X-ray peaks are aligned with the E > 100 MeV data 
within Acf) ^ 400 /is (a value obtained fitting the peaks 
with Gaussians) providing an additional test of the AG- 
ILE phasing accuracy. 

The effective time resolution of AGILE light-curves 
results from the combination of the different steps in- 
volved in the processing of gainma-ray photon arrival 
times. The on-board time tagging accuracy is a mere 
^1 /is, with negligible dead time. For comparison, the 
corresponding EGRET time tagging accuracy was ~100 
/Its. The precise GPS space-time positioning of AGILE 
spacecraft (Argan et al. 2004) allows for the transforma- 
tion from UTC to Solar System barycenter time-frame 
(TDB) with only a moderate loss (<10 fis) of the intrinsic 
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Fig. 8.— PSR B 1706-44 light-curves (P ~ 102.5 ms) for different 
energy bands {E > 30 MeV, G class events, 40 bins, resolution: 
~2.6 ms; 0.1 < i? < 1 GeV, G+L class events, 20 bins, resolution: 
~5.1 ms; E > 1 GeV, G+L class events, 10 bins, resolution: ~10.2 
ms) obtained by integrating all available data presented in Table 1. 
The radio ephcmcris and the 1024 bins light-curve (bottom panel) 
are obtained by the analysis of 20 ToAs observed by Parkes radio 
telescope (observation interval 54,220-54,562 MJD). 



instrumental time accuracy. The innovative folding tech- 
nique described in § 4, accotmting for pulsar timing noise, 
also reduces smearing effects in the light-curves, fully ex- 
ploiting all the information from contemporary radio ob- 
servations. In summary, the effective time resolution of 
the current AGILE pulsar light-curves (and then multi- 
wavelength phasing accuracy assuming terr =^ 1 fJ-s) is 
mainly limited by the available count statistics and can 
be estimated by: 



P _ a'{C^ + 2B) 



(4) 



where P is pulsar period, TV is the number of bins in 
the light-curve histogram, a is the signal-to-noise ra- 
tio, Cp are the pulsed counts and B are the background 
counts. In order to keep the average signal-to-noise ra- 
tio of light-curve bins (during the on-pulse phase) at a 
reasonable level (>3cr), the resulting effective time res- 
olution is constrained to 200 500 /is. At present, the 
best effective time resolution (^200 fis) is obtained for 
the 400-bin light-curve of Vela (G+L class selection) al- 
though a 100-bins light-curve (Figure 5) is better suited 
to study pulse shapes and to search for possible weak 
features. The effective time resolution will obviously im- 
prove with exposure time At oc and a resolution ;$50 



/xs is expected after two years of AGILE observations of 
Vela. 

6. DISCUSSION 

With about 10 years since the last gamma-ray obser- 
vations of Crab, Vela, Ceminga (Fierro et al. 1998) and 
B 1706-44 (Thompson et al. 1996) by CGRO, the im- 
proved time resolution of AGILE and the much longer 
observation campaigns in progress are now offering the 
possibility of both to search for new features in the shape 
of the light-curves of these gamma-ray pulsars and to 
investigate the possible occurrence of variations in the 
gamma-ray pulsed flux parameters. 

After nine months of observations in the frame of the 
Science Verification Phase (2007 July-November) and 
of the Scientific Pointing Program AO 1 (pointings 1- 
7, 2007 December 1-2008 April 10), AGILE reached an 
exposure (E>100 MeV) of the Vela region >10^ cm^ s 
(^10,000 pulsed cotmts from Vela), comparable to that 
of the nine-year life of EGRET (although AGILE data 
have a higher residual particle background), and an even 
better exposure (1.5 x lO'' cm^ s) in the core region of 
the Galactic Plane {I = 310° 340°) corresponding to the 
Southern Hemisphere . In fact, AGILE observed ~2,400 
pulsed counts from PSR B1706— 44 up to date, a factor 
1.5 better counts statistics than EGRET for this source. 
For Crab and Ceminga, an exposure level comparable 
with that obtained by EGRET will be reached at the 
end of the AO 1 pointing number 15 in the Anti-Center 
region (October, 2008). 

6.1. The gamma-ray light- curves 

The plots shown in Figures 5, 6, 7, and 8 allow us to 
start assessing new features in gamma-ray pulsar light- 
curves. Narrower and better resolved main peaks are re- 
vealed, together with previously unknown secondary fea- 
tures, to be confirmed when more count statistics (and an 
improved particle background filtering) will be available. 

The Vela light-curves for different energy bands are 
shown in Figure 5. A Gaussian fit to the Vela main peak 
(PI) atE> 100 MeV provides a FWHM of 0.018 ±0.002 
centered at = 0.1339 ± 0.0007 in phase consistent with 
the EGRET observations (Kanbach et al., 1994). The 
apex of the main peak is resolved by AGILE with a width 
of -^0.8 ms and its apparent trail (^ = 0.13-0.3) could be 
due to the occurrence of one or more secondary peaks. In 
fact, marginal evidence of a relatively narrow lower peak 
(P3) at ~ 0.25 is present in the E <100 MeV (3.3cr 
fluctuation with respect to the average interpulsc rate) 
and in the E >1 GeV (3.8cr) light-curves. P3 is located 
at the phase of the optical peak 1, and also at the phase 
of a bump predicted in a two-pole caustic model (Dyks 
et al. 2004), due to overlapping field lines from opposite 
poles near the light cylinder. In the outer gap model, this 
bump is the first peak in the light-curve, and also comes 
from very near the light cylinder. The peak at = 0.5- 
0.6 (P2) in the E > 100 MeV hght-curve cannot be satis- 
factorily fit with a single Gaussian or a Lorentzian curve 
{Xr > 3), due the possible presence of a bump at (/) ^ 0.5 
(P2a). A fit with two Gaussian provides 4> = 0.560±0.001 
(FWHM 0.031 ± 0.003) and (f> = 0.49 ± 0.01 for the 
phases of the major peak (P2b) and lower peak apex 
(P2a), respectively. The phase separation between the 
main gamma-ray peaks A(p = 0.426 ± 0.002, as well as 
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that between the gamma-ray and the radio peak, arc 
unchanged since EGRET observations (Kanbach ct al. 
1994; Ramanamurthy et al. 1995). In the E < 100 MeV 
band a secondary peak+valley structure (P4) appears 
at ^ ~ 0.9 (~4iT level). It is worth noting its symmet- 
ric position around the radio peak with respect to main 
peak (PI) and a possible correlation of P4 with features 
seen in the X-ray light-curves (Manzali et al. 2007); P4 
also coincides with peak number 3 of the RossiXTE light 
curve (Harding et al. 2002). 

Geminga shows an £^ > 100 MeV light-curve (Figure 6) 
with properties similar to those of Vela. Apart from the 
major peaks P2 (</> = 0.999±0.002, FWHM 0.062±0.008) 
and PI (with apex = 0.507±0.004, FWHM 0.08±0.01), 
secondary peaks are seen trailing PI and leading P2: the 
main peak's P2 leading trail {(j) = 0.8-1) could be pos- 
sibly associated to unresolved multiple peaks while PI 
displays a "bump" at ~ 0.55-0.6. PI seems to be char- 
acterized by a double structure: a fit with a simple Gaus- 
sian yields Xr 1-2 to be compared with a double Gaus- 
sian model providing < 1. This feature seems still 
present even considering different observation block sep- 
arately. The main-peak separation {A(j> = 0.508 ± 0.007 
at E > 100 MeV) is greater than the value obtained for 
Vela and it decreases slightly with energy. The 2-10 keV 
X-ray light-curve shows a peak (P3) in correspondence 
with a possible excess in the hard gamma-ray band and 
a broad top-hat shaped feature partially overlapping in 
phase with PI. 

The Crab hght-curve for E > 100 MeV (Figure 7) 
shows a previously unknown broad feature at (/> = 0.65 — 
0.8 (P3), in addition to the main peaks PI and P2 {(f)pi = 
0.999±0.002, FWHM 0.054±0.005: cf>P2 = 0.382±0.008, 
FWHM 0.14 ± 0.04). The probabihty of P3 being a 
background fluctuation is of ~10~** (3.7cr). P2 could 
be possibly resolved in two sub-peaks in future longer 
observations. P3 is coincident with the feature HFC2 
that appears in the radio profile above 4 GHz (Moffett 
& Hankins 1996). From the polarization of this compo- 
nent, Moffett & Hankins (1999) suggest that this peak 
may come from a lower emission altitude, near the polar 
cap. P3 is actually at a phase that could plausibly come 
from low altitude cascades in the slot gap model (Mus- 
limov & Harding 2004, Figure 2), the pairs from which 
may be also causing the HFC2 radio component (s), while 
PI and P2 comes from the high-altitude slot gap. 

According to the observations of SAS-2, COS B and 
CGRO/EGRET, the ratio P2/P1 of the main peak in- 
tensities could present a variability pattern (possibly as- 
cribed to the nutation of the neutron star) that can be fit 
with a sinusoid with a period of ~13.5 years (Kanbach 
1990; Ramanamurthy et al. 1995) although this is not 
required by EGRET data alone (Tompkins et al. 1997). 
We observed a P2/P1 intensity ratio 0.66 ± 0.10 in good 
agreement with the value of ~0.59 predicted for 54,350 
MJD (for the energy range 50 MeV-3 GeV). Unfortu- 
nately, our P2/P1 value is similar - within the errors - to 
the EGRET determination (~0.5): then an unambigu- 
ous assessment of the origin of this possible phenomenol- 
ogy will require measurements close to the epoch (56,150 
MJD) corresponding to the predicted maximum or the 
intensity ratio P2/P1 (~1.4). Variability should also be 
invoked to explain the possible detection of P3, which 
was never seen before in the EGRET database in spite of 



an overall exposure comparable to that reached by AG- 
ILE so far. We note that the main peak intensity ratios 
computed for Vela (P2/P1 = 0.91 ± 0.07) and Geminga 
(P2/P1 = 0.8 ± 0.1) do not yield evidence of significant 
variations with respect to past observations (Ramana- 
murthy et al. 1995). 

The AGILE light-curves of PSR B 1706-44 are shown 
in Figure 8. The broad-band hght-curve {E > 30 MeV) 
clearly shows two peaks {(j)pi = 0.211 ± 0.007, (/)p2 = 
0.448 ± 0.005) bracketing considerable bridge emission 
(contributing to >50% of the pulsed counts) while in the 
0.1 1 GcV range the peaks cannot be discerned from the 
bridge emission and the pulsar profile presents an unre- 
solved broad (A0 = 0.3 - 0.4) single peak. PSR B 1706- 
44 is a young (~ 2 x 10* yr) and energetic (3.4 x 10^^ erg 
s~^) 102.5 ms pulsar (Johnston et al. 1992) with emis- 
sion properties similar to Vela (Becker & Pavlov 2002). 
Double peaked PSR B 1706-44 is then in fact "Vela-like" 
not only energetically, but also with respect to the off- 
set between the maxima of the high energy and the radio 
profiles, with neither of the two gamma-ray narrow pulses 
aligned to the radio peak. 

AGILE allowed for a long monitoring of gamma-ray 
pulsar light-curves shapes. We carefully looked for pos- 
sible lightcurve variations by KS tests (two-dimensional 
Kolmogorof-Smirnoff test [KS]; Peacock 1983; Fasano & 
Franceschini 1987). For each pulsar, different gamma- 
ray light-curves (with 10, 20 and 40 bins) were ob- 
tained grouping contiguous data set and requiring at 
least 30 counts bin~^ (300-1000 counts for each light- 
curve). Each light-curve was compared by KS test with 
each other and with the average shape corresponding 
to the entire data set. No pulse shape variation was 
detected with a significance >3 cr on timescales ranging 
from 1 days (Vela) to few months (Crab, Geminga and 
PSR B 1706-44). 

6.2. Implications for the emission models 

Pulsars derive their emitting power from rotational en- 
ergy loss owing to the relativistic acceleration (up to 
Lorentz factor F ~ 10^-10^) of charged particles by very 
high electric potentials induced by the rotating magnetic 
fields. The charge density that builds up in a neutron 
star magnetosphere is able to short out the electric field 
parallel to the magnetic field everywhere except a few 
locations of non-force-free "gaps" . It is unclear whether 
these acceleration regions can form in the strong field 
near (<1 stellar radii) the neutron star surface (polar 
cap/low-altitude slot gap model: Daugherty & Harding 
(1996) and Muslimov & Harding (2003); high-altitude 
slot gap model: Muslimov & Harding (2004) and Harding 
et al. (2008)), in the outer magnetosphere near the speed 
of light cyhnder (outer gap model: Cheng et al. (1986); 
Romani (1996); Hirotani & Shibata (2001); Takata & 
Chang (2007)) or even beyond in the wind zone (Petri & 
Kirk 2005). 

In polar cap/low altitude slot gap models, gamma-rays 
result from magnetic pair cascades induced by curvature 
or inverse-Compton photons. The spectrum is domi- 
nated by synchrotron radiation of pairs at lower energies 
(<1 GeV) and by curvature radiation at higher energies 
(>1 GeV). In the high-altitude slot gap, gamma-rays re- 
sult from curvature radiation and synchrotron radiation 
of primary electrons (no cascade). In outer gap models. 
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gamma-rays result from photon-photon pair cascades in- 
duced by curvature radiation. The spectrum is domi- 
nated by pair synchrotron radiation below 20 MeV, cur- 
vature radiation above 20 MeV and inverse Compton at 
1-10 TeV. 

The size and spectrum of the emitting regions are then 
directly related to the intensity and location in the mag- 
netosphere of the accelerating electric fields. Strong fields 
imply thin accelerator gaps, while weaker fields are as- 
sociated to thick gaps: the acceleration zone grows big- 
ger as the particle must accelerate over larger distance 
to radiate pair-production photons. On an observational 
point of view, it is in turn expected that acceleration gap 
sizes are related to the width of light-curve peaks. The 
highly relativistic particles emit photons at very small 
(~l/r) angles to the open magnetic field lines. The the- 
oretical width of a light-curve peak associated to an in- 
finitely small gap would be then At = P/(27rr), a value 
typically smaller than 1 iis. Therefore, the width of the 
apex of the peak can be related to the core gap size. 
For example, the At ~ 800 /xs width of the Vela pulsar 
peak (PI) resolved by AGILE implies a projected gap 
core width of ~1 km for gap height of ~1 stellar radii. 
Broader peaks as P2 and P3 in the Crab light-curve in- 
volves instead magnetospheric region tens of kilometers 
long. The relation of the width of the peaks to the ac- 
celeration gap size works for the low altitude emission 
models but is complicated by relativistic effects (aberra- 
tion, retardation) in the high-altitude emission models, 
where the peaks are formed by caustics (Dyks & Rudak 
2003). The peak widths will still depend on the width of 
the acceleration (outer or slot) gap, but one must per- 
form a comparison with detailed models to constrain the 
physical size of the gap. 

The likely presence of multiple contiguous peaks as 
P2a-b in the Vela light-curve, PI in Geminga and pos- 
sibly within the wide P2 broad peak in Crab could 
be related to short-term oscillation of gap locations. 
Light-curve variations on timescales >1 days are any- 
way excluded by the KS tests described above. Alter- 
natively, the apparent superimposition of different gaps 
placed at different heights in the magnetosphere could 
be another plausible explanation for multiple-contiguous 
peaks, while the hypothesis of jumps in the acceleration 
electric fields strength would lack of evident physical jus- 
tification. In general, the presence of multiple (contigu- 
ous or not) structures in the light-curves is difficult to 
fit in a scenario alternatively involving polar cap, outer 
gap or wind zone models exclusively. The particle ac- 
celeration may well be simultaneously occurring in all 
the regions predicted by these models. Indeed, a three 
peak light-curve can still be explained invoking polar cap 
gaps alone provided that both polar caps can cross the 
line-of-sight. 

Wherever in the magnetosphere a gap can form, a 
multi-frequency emission may occur along an hollow cone 
(due to magnetosphere symmetries) or any other suitable 
surface: light-curve peaks are generated when the view- 
ing angle from any given location on this surface to the 
observer also crosses the gap. The pulse spectrum slope 
depends on the accelerating field as well as on the mag- 
netic fields, which strongly affect the synchrotron emis- 
sion efiiciency and the pair attenuation along the cas- 
cade path: the stronger the magnetic field is, the steeper 



the spectrum is. An outer gap can in principle gener- 
ate pulse spectra extending up to tens of GeV (Zhang & 
Cheng 2000), while gaps close to surface high magnetic 
fields can produce steep spectra (photon index <3; Ra- 
manamurthy et al. 1996) and cut-off in the tens of MeV 
range (PSR 1509-58, Harding et al. 1997). The absence 
of corresponding gamma-ray pulse in phase with radio 
main peaks in Vela and PSR B 1706 -44 could be then 
ascribed to gap-dependent spectral slope and not only to 
beam angle and viewing geometry difi^erences. 

The multiplicity and variety of features seen in AG- 
ILE light-curves can pave the way to a parameterized 
standard model (e.g. with adjustable accelerating elec- 
tric fields strength, and location in the magnetosphere) 
for pulsars gaps and their corresponding observed high- 
energy pulses. In this perspective, the AGILE light- 
curves time resolution, currently limited only by the 
(continuously increasing) source counts statistics, will 
eventually yield a pulsar gaps map by coupling timing 
analysis and phase resolved spectral analysis. 

6.3. The Vela glitch of August 2007 

During early AGILE observations, Vela experienced a 
weak glitch clearly detected in radio as a discontinuity in 
the pulsar's spin parameters. 

Glitches are small {Avjv ~ 10~^-10~^) and sudden 
(<1 day) discontinuous increases in the pulsar frequency, 
often followed by a recovery (1-100 days) to the pre- 
glitch frequency. About ~6% of pulsars are known to 
have shown glitches^^), with a higher incidence of events 
in younger pulsars. This phenomenon is potentially a 
very promising tool for probing the physics of the neu- 
tron star interiors (Lyne et al. 2000). Although no gen- 
eral consensus has been reached to date about the origin 
of the glitches, many models are based on the exchange of 
angular momentum between the superfluid neutron star 
core and its normal, solid crust (Ruderman 1976, 1991; 
Alpar et al. 1984b, a). This angular momentum trans- 
fer may excite starquake waves, propagating toward the 
neutron star surface. Since the magnetic field frozen in 
the crust is "shaken", the resulting oscillating electro- 
magnetic potential could generate strong electric fields 
parallel to the magnetic field, which in turn would accel- 
erate particles to relativistic energies, possibly emitting 
a burst of high-energy radiation. 

Since the first observation of a pulsar glitch in 1969 
(Radhakrishnan & Manchester 1969), Vela has shown 
^10 major glitches. Due to its large field of view, the 
quest for possible gamma-ray bursting behavior due to 
a glitch is then an effective opportunity for AGILE. De- 
spite the fact that the August 2007 glitch is a weak one 
[Aviv ^ 10"''), it is worthwhile to search for a signal in 
the AGILE data. 

The characteristic energy of a pulsar glitch can be 
roughly estimated from the associated pulsar frequency 
jump E'giitch = AErot = 47r^7i/Az/, where I (~10^^ g 
cm^) is the neutron star momentum of inertia. The cor- 
responding expected gamma-ray counts would be: 



-E'glitch^eff 

And'^E^ 



(5) 



where r] = [0, 1] is the unknown conversion efiiciency of 
See e.g. http://www.atnf.csiro.au/research/pulsar/psrcat/. 
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Fig. 9. — Unfolded Vela pulsar light-curve (4.5 min binsize, 
i? > 50 MeV, G-class events). Dashed lines bracket glitch epoch 
uncertainty range (54, 313 ± 3 MJD). A >5<T count excess at 
~54,312.693 MJD could be associated to gamma-ray bursting emis- 
sion from the glitch. 

the glitch energy to gamma-ray emission, d (^0.3 kpc) 
is the pulsar distance, E-^ (^300 McV) is the average 
gamma-ray photon energy assuming a spectral photon 
index F = — 2, and Aoff is the AGILE effective area. 

Even in the virtual limit assumption that the entire 
glitch energy could be driven into gamma-ray emission, 
a weak glitch with a frequency shift of Au/u = 1.3 x 10~^, 
as that observed in August 2007, cannot produce a strong 
signal (Cfi^^h ^ 100-200 counts), if the core fluence is 
spread in ~1 day. In fact, no excess on daily timescales 
was detected, although for much shorter timescales of 
3 6 minutes a >5<t excess (^15 counts) in the photon 
counts is seen at -54,312.693 MJD (Figure 9). 

On the other hand, stronger Vela glitches, as that of 
1988 Christmas (McCulloch et al. 1990) with a frequency 
shift of Av/v = 2 X 10~^, could in principle produce more 
significant transient gamma-ray emission. Typical count- 
rate from Vela is '--^100-200 counts day^^, then a fluence 
of >1000 counts in ~l-day or less from the hypothetical 
gamma-ray glitch burst should be easily detectable. Ac- 
cording to Equation (5), such a flux could arise from a 
glitch with Avjv ^ 10~^ (typical Vela glitch size), con- 
verting a relatively small fraction {r\ ~ 0.1) of its energy 
in gamma-rays. The chance occurrence of a strong Vela 
glitch in the A GILE field of view over three years of mis- 
sion operations is of ~20%. 

7. CONCLUSIONS 

AGILE collected ~15,000 pulsed counts from known 

gamma-ray pulsars during its first 9 months of opera- 
tions. The AGILE PDHU clock stability, coupled with 
the exploitation of pulsar timing noise information, al- 
lows for pulsar period fitting with discrepancies with re- 
spect to radio measurements at the level of the period 
search resolution (~10~"'^^ s) over the long gamma-ray 
data span (>6 months). Thanks to AGILE GPS-based 
high time tagging accuracy ('^l /is), the effective time 
resolution of AGILE light-curves is limited only by the 
count statistics at current level of exposure (~1.2 x 10® 



cm^ s for the Vela region). The best effective time reso- 
lution obtained for Vela observations is of ^200 /is for a 
signal-to-noise >3 a in the on-pulse light-curve bins. An 
improved effective time resolution < 50/US is expected af- 
ter three years of AGILE observations. 

AGILE multi- wavelength phasing of the four gamma- 
ray pulsars is consistent with the results obtained by 
EGRET, although the high resolution AGILE hght- 
curves show narrower and structured peaks and new in- 
teresting features. In particular, a third peak is possibly 
detected at ^3. 7 a level in the Crab light-curve and sev- 
eral interesting features seem present in the Vela light- 
curves. 

In any case, the highly structured light-curves hint at 
a complex scenario for the sites of particle acceleration 
in the pulsar magnetospheres, implying different electric 
gaps with physical properties probably mostly related to 
their height above the neutron star surface. Alterna- 
tively, slight spatial oscillations of the gap locations on 
timescales <1 day could be invoked to explain the mul- 
tiple contiguous peaks seen in the light-curves. 

The foreseen balance of count statistics at the end 
of AO 1 observing program will allow for phase-resolved 
spectral analysis of the light-curves and, correspondingly, 
a spatial mapping of the magnetospheric gaps (their al- 
titude above the neutron star surface being possibly re- 
lated to the shape of their spectra). 

We finally note that the timing calibration and tests 
presented in this paper pave the way to an effective search 
for new gamma-ray pulsars with AGILE. The negligible 
discrepancies among the radio and the gamma-ray pul- 
sar spin parameters seen in the known gamma-ray pul- 
sars imply that the direct folding of the AGILE data 
on new pulsar candidates is a safe and efficient proce- 
dure, when the folding parameters are obtained from 
radio/X-rays ephemeris having suitable epoch range of 
validity. In fact, gamma-ray period search trials around 
radio/X-rays solutions will be unnecessary when simul- 
taneous ephemeris are available, strongly improving the 
detection significance for faint sources. According to the 
predicted gamma-ray pulsar luminosities {L^ oc E^^'^/(P, 
e.g. PeUizzoni et al. 2004), the detection with AGILE of 
top-ranked Vela-like pulsars is then expected as soon as 
exposure levels >10'' cm^ s (£■ > 100 MeV) will be at- 
tained for clean G-class events. 
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